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Abstract 

For a high-quality single crystal of UPd3 we present the relevant elastic 
constants and ultrasonic attenuation data. In addition to the magnetic phase 
transition at T2 = 4.4±0.1K and the quadrupolar transition at T\ ~ 6.8K, we 
find orbital ordering at To = 7.6 =L O.IK concomitant with a symmetry change 
from hexagonal to orthorhombic. A striking feature is the splitting of the 
phase transition at T\ into a second-order transition at T + \ = 6.9±0.05K and 
a first-order transition at T_i = 6.7 ± 0.05K. For the four phase transitions, 
the quadrupolar order parameters and the respective symmetry changes are 
specified. 

PACS Numbers: 71.45.-d, 62.20.Dc, 71.70.Ms 

The actinide intermetallic nonheavy fermion system UPd3 and the heavy fermion super- 
conductor UPt3 are hexagonal compounds with space group PQ^/mmc, exhibiting magnetic 
phase transitions at T 2 = 4.4K and ~ 5K, respectively. This finding is very surprising, 
because in both compounds the ordered magnetic moment is extremely small (/x ~ lO -2 //^). 



1 



Furthermore, the 5/ electrons of UPt3 are found to be predominantly itinerant, whereas 
UPd 3 belongs to the rather few metallic materials where well-localized 5/ electrons and 
ionic long-range quadrupolar ordering at T\ = 6.8K has been identified. The possibility 
that quadrupolar order occurs at T\ was first recognized in specific heatJi and thermal ex- 
pansion! measurements, and the strong quadrupole-quadrupole interaction between the U 
ions was confirmed for the first time in sound velocity experiment M, where also the rel- 
evant quadrupole components contributing to the magnetoelastic coupling were specified. 
The ultimate proof of the quadrupolar nature of the phase transition, however, was given 
in neutron scattering experimental, where the long-range lattice distortions which always 
accompany§ a quadrupolar phase transition were determined. 

The main reason for the interest in the quadrupolar phase transition of UPd3 is that the 
knowledge of the ordering components of the quadrupole tensor not only provides detailed 
information about the long-range ordering of the 5/ 2 electronic charge distribution of the 
uranium U 4+ ions (ground state 3 H 4 ), but also of long-range spatial correlations of the 
orbital degrees of freedom. This is so, because the average values of the components of 
the electric quadrupole tensor Q are not only a measure of deviations of the localized 5/ 
uranium electrons from a spherical charge distribution but, according to the identity oc 
{3(JjJj + JjJi)/2 + J(J + 1)}, are also a measure of correlations between the components 
Jj and Jj of the total angular momentum J. We mention that long-range quadrupolar (i.e., 
orbital) ordering is also of current interest in the context of the physics behind the colossal 
magnetoresistance in the manganatesi or the spin-Peierls and metal-insulator transitions!! 
in the vanadates. 

In order to deduce from the lattice distortions the ordering quadrupole components Qij 
(i.e., the order parameters), the complete strain tensor e or, more strictly speaking, the 
symmetry adapted strains Ey must be known because the latter are proportionali to the 
average value of the symmetry adapted components of the quadrupole tensor. In neutron 
diffraction experiments, however, the full strain tensor is difficult to deduce from the available 
Bragg peaks. This is one of the reasons why details concerning the various phase transitions 



and the proper order parameters in UPd3 are still under discussion. Also, it is still an open 
question to what extent the inter-ion quadrupole-quadrupole interaction is responsible for 
the three reported phase transitions found at T2 = 4.4K, T\ = 6.8K and To = 7. 8KB. 

It is well known that ultrasound is well suitedi to answer these questions because the 
coupling of the acoustic strain e to the ionic quadrupole moment tensor Q is, in general, large 
and the coefficients of the elastic stiffness tensor C depend in a characteristic manner!'! (see 
Table |) on the components of the so-called quadrupolar strain susceptibility x ■ The latter 
is defined as SQ = x V and therefore is a measure of the response 8Q of the quadrupole 
moment tensor Q to the strain-induced electric field gradient tensor V at the uranium sites. 
The components are related to the elastic strain tensor e via = J2 Sijki £ kh where 
the fourth-rank tensor S is the so-called field gradient elastic strain tensor which may be 
determined! in nuclear acoustic resonance experiments: this should not be confused with 
the elastic compliance tensor. The quadrupolar contribution Cp to a symmetry elastic 



constant Cp is given b; 
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where Ny is the number density of quadrupoles, g' T the two ion quadrupole-quadrupole 
coupling constant and (depending on the irreducible representation T), is a linear 
combination of the S-tensor components for g' T = (i.e., in the absence of the inter-ion 
quadrupole-quadrupole interaction). In terms of fluctuations [i.e., of the variance (AQ r ) 2 } 
of the symmetry adapted quadrupole tensor components Qr, the quadrupolar strain suscep- 
tibility Xr may be written in the thermostatic (zero-frequency) limit as Xr ^ (^Cjr) 2 / 
(fcgT), where, in general, five relevant quadrupole tensor components exist which, in the 
symmetry adapted form, we denote by Q zz , Q x 2_ y 2 = Q xx - Q yy , Q xy , Q yz , and Q zx . If 
Qr is the order parameter, then a pronounced decrease of the related elastic constants Cp 
(see Table |) should be observed at a second-order quadrupolar phase transition, because 
in that case, the fluctuations of Qr (and accordingly the strain susceptibility Xr ) should 
become large. Provided that at a quadrupolar phase transition the formation of domains or 



its influence on the elastic behavior can be minimized such as to play only a minor role, then 
it becomes immediately evident that ultra-sound should be well suited to study long-range 
quadrupolar (i.e., orbital) ordering. 

The efficiency of ultrasound is demonstrated in Figs. 1(a) and 2(a) for longitudinal acous- 
tic waves propagating along the c axis of a rectangular single crystal of UPd3 (with the linear 
dimensions 3.6 x 5.2 x 7.3 mm 3 ). As can be clearly seen, both the attenuation (phonon loss 
rate) and the sound velocity of the C33 mode display three distinct phase transitions at 
T 2 = 4.4 ± O.IK, Ti « 6.8K and T = 7.6 ± O.IK, whereas the specific heat0 and magnetic 
susceptibility^ of samples of the same batch exhibit pronounced peaks only at T\, and T\ 
and T 2 , respectively. Previous sound velocityi'El and neutron scattering experiments^ re- 
vealed antiferroquadrupolar ordering at T\ (Refs. 13 and 14) on the quasi-cubic uranium 
sites and it was thought that the concomitant structural transformation is from hexagonal 
to trigonal00 or to orthorhombid or monoclinici. At first sight, it is surprising that this 
transition dominates the specific heat but (compared to the magnitude of the ultrasound 
absorption peaks observed at T 2 and T ) contributes only little to the phonon-loss rate. Ac- 
cording to Table | and to Refs. 5 and 15, on the other hand, (Q zz ) is not the proper order 
parameter of a quadrupolar phase transition concomitant with a structural transformation 
from hexagonal to trigonal (monoclinic or orthorhombic) . We therefore do not expect that 
the phonon-loss rate of the C33 mode will be modified essentially in the vicinity of T\ which 
is in agreement with our experimental findings (see Fig. 1). The absence of a huge absorp- 
tion peak at Ti (see Fig. 1) of both the C33 and the C44 mode further indicates that at 
T\ the contribution of incoherent strains (originating from domains or domain walls) may 
be considered as negligibly small. This suggests that the ordered phase is the quadrupolar 
triple-q stated, because the triple-q state avoids the formation of stochastically distributed 
domains. The development of the triple-q state is associated with modulations of the elec- 
tronic charge distribution. These modulations are a superposition of three plane waves with 
wave vectors qi, q 2 and q3, where the wave vectors q 2 , and q3 are obtained from q x by 
rotations of ±27r/3. The trigonal triple-q phase (space group P3ml) was discovered^ in 



neutron diffraction experiments and, within experimental uncertainty, the respective phase 
transition at Ti was found to be continuous^. As we will see below, however, this is in 
contradiction to our findings. 

We have mentioned already that the average value (Q zz ) cannot serve as a primary 
order parameter because it is already nonzero above Ti (and T ). The formation of the 
triple-q phase may, nevertheless, modify the C33 mode because the triple-q state alters^ 
the Q zz moments of the cubic site ions. Surprisingly (see Fig. 2) the C33 mode does not 
soften at Ti upon cooling (as expected for a continuous phase transition), but shows a 
steplike increase which suggests that the transition at Ti is most likely of first order. Even 
more surprising is the finding [see inset of Fig. 1(a)] that at around Ti not one but two 
clearly distinguishable absorption peaks (and accordingly two phase transitions) appear at 
T_i = 6.7K and T + \ = 6.9K, where the stiffness coefficients C33 and C44 turn out to be 
hysteretic in the vicinity of T_i (see Fig. 3) but not at T+\. In addition, [see Fig. 1(b)] 
we also have observed at T_i a pronounced hysteresis in the attenuation of the C44 mode. 
The quadrupolar transition into the trigonal triple-q state must therefore be of first order, 
whereas the transition at T +1 is most likely of second order. We would like to emphasize 
that, to the best of our knowledge, for UPd3 neither ultrasound absorption measurements nor 
the existence of the two neighboring phase transitions at Ti have been previously reported. 
Compared to other experimental techniques, it is furthermore worthwhile to note (see Figs. 1 
and 2) that at zero magnetic field the phase transition at T = 7.6K is most apparent in 
ultrasound. 

So far, we have confined ourselves mainly to the C33 mode which only reflects the tem- 
perature dependence of the secondary order parameter (Q zz ). For hexagonal symmetry (see 
Table |) a softening of the C 44 mode is attributed to an ordering of the quadrupole com- 
ponents (Q yz ) and/or (Q zx ), whereas the order parameters related to C 66 = (Cn — Ci 2 )/2 
(or to C11 and C22) are the degenerate quadrupoles (Q xy ) (or (Q x a -y a ))- We now consider 
the temperature dependence of C44 and Cqq. As can be seen in Fig. 2, with decreasing 
temperature, C44 exhibits a small dip at T_i, followed by a steep decrease at the magnetic 
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transition temperature T2 = 4.4K which is accompanied by a very sharp peak (see Fig. 1) 
in the phonon loss rate. However, the most pronounced softening of the Cqq mode appears 
at To = 7.6K, not at T\ as reported for a sample of much lower crystal quality and size by 
other authorsi. The softening of the C44 and modes, which starts already far above the 
transition temperatures T 2 and T , is very strong and it has been shownlll that below 150K 
the temperature dependence of the elastic constants is essentially due to the crystal-field 
splitting at the quasicubic uranium sites. We therefore conclude that not only at Ti, but 
also at To and at the magnetic transition temperature T2, the driving force of the respective 
phase transitions is of quadrupolar origin. Hence, the quadrupolar order parameter of the 
magnetic phase (T < T 2 ) is most likely (see below) a linear combinational of (Q yx ) and 
(Qx 2 - y 2 ), whilst at T the relevant order parameters are the degenerate quadrupoles (Q xy ) 
and (Qx 2 - y 2 )- Here, degenerate means that the response of the respective quadrupoles is 
specified by the same strain susceptibility. 

The degenerate order parameters (Q x 2_ y 2) and (Q xy ) which, according to Table |, will 
modify Cu, C22 and Cqq = (Cu — C12V2, are associated with a transition from hexagonal 
to orthorhombic0'0 or to monoclinic0, whereas the order parameters (Q yz ) and (Q zx ) will 
change C44 and C55 and entail a transition from a hexagonal to either a monoclinic or triclinic 
phaseSB. Here, the formation of a monoclinic or triclinic phase depends on whether or not 
both order parameters are present simultaneously. Since no anomaly is seen at T = 7.6K 
in the elastic modulus C44, whereas the Cu and C22 modes (which are not shown here) and 
the C*66 mode soften significantly, it follows that at T we pass with decreasing temperature 
from the hexagonal to an orthorhombic or monoclinic phase. When analyzing the phase 
transitions below To, we therefore have to take into account that for T + i < T < To, the 
symmetry is orthorhombic or monoclinic but not hexagonal (which in the literature was 
taken for granted until recently). The monoclinic phase, however, may be excluded because 
in that case the only possible second-order transition at T +1 would hM from a monoclinic to 
a triclinic phase, which is in contrast to our experimental findings (see Fig. 3), because such 
a transition should alter C44 significantly. At Tq the symmetry of UPd3, therefore, passes at 



decreasing temperature from hexagonal to orthorhombic and the only possible-second order 
transitionliUll at T +1 is from the orthorhombic to a monoclinic state. At T +1 we therefore 
have to consider0@ (Q yz ), (Qzx), and (Q xy ) as possible order parameters, where (Q yz ) and 
(Qzx) ma y be excludedl^E because C44 (see Fig. 3) and C55 do not alter significantly at T + \. 
Concerning (Q xy ), however, which should modify0 Cqq, we cannot prove experimentally 
whether or not (Q xy ) becomes an order parameter at T+i, because in the temperature range 
between 6.6K and 7.8K, the attenuation of the Cqq mode becomes too large. We note that 
in agreement with Ref.0 and our experimental findings, the phase transition at TLi (i.e., 
from the monoclinic to the trigonal triple-q state) is of first order. The softening of the C44 
mode at T 2 = 4.4K further indicates^! that at decreasing temperature, UPd3 transforms at 
T 2 from the trigonal triple-q state into a monoclinic (magnetic) state, where the transition is 
expected^ to be of first order. Within experimental resolution, however, the C44 mode does 
not show a signature of a first-order transition which indicates that at T 2 it is possibly not 
enough to consider the magnetoelastic coupling only. According to Ref. 15 the transition 
from the triple-q into the monoclinic phase should be accompanied by the loss of the ternary 
axis and the order parameter should be a linear combination of (Q yz ) and (Q x 2_ y 2). It is 
there-fore expected that both C44 and Cqq will soften at T 2 , which is confirmed by our 
experiments (see Fig. 2). In the immediate vicinity of T 2 , however, the attenuation of the 
C*66 mode becomes too large to enable a more detailed analysis of this phase transition. 

Based upon the crystal-field split (J = 4) multiplets determined by Buyers et a/.,0 
the temperature dependence of C 44 has been evaluated previously in Ref. 13 yielding the 
quadrupolar coupling constants Nv(S^) 2 /Cj$ = 0.6K and g' 44 = — 6K. Since the ordered 
magnetic moments in UPd3 are so extremely smallB and the quadrupolar coupling constants 
are so large, the pronounced softening of the C 44 mode indicates that the quadrupole- 
quadrupole interaction is the most probable driving mechanism of the magnetic phase tran- 
sition at T 2 . Owing to the negative sign of the quadrupole coupling constant g 44 we further 
expect0@ antiferroquadrupolar ordering at T 2 and the formation of a very complex an- 
tiferromagnetic state, which is currently under investigation and will be the subject of a 



forthcoming paper. Summarizing our results, we have shown that ultrasound is a powerful 
tool for investigating orbital ordering and that the large variety of interesting phase transi- 
tions in UPd3 is of quadrupolar origin. We hope that our findings may also offer some new 
ideas in the attempt of a better understanding of the formation of the heavy fermion state 
and the antiferromagnetic transition in UPt 3 . 

This work was supported by the Deutsche Forschungsgemeinschaft (Grant No. MU 696/ 
2-2). We acknowledge many helpful and clarifying discussions with K.-H. Hock. 
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FIG. 3. Hysteretic behavior of the C33 and C44 elastic mode at about 22MHz. 
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TABLES 

TABLE I. Quadrupolar contribution C^g to the elastic constants in Voigt notation, the 
strain-induced quadrupole interaction Hamiltonian Hq and the relevant strain susceptibilities Xr"^ 
for hexagonal crystal symmetry and different acoustic strain modes £jj. The indices x, y, and z 
refer to the a-, b-, and c-axis of the crystal. 
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